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We  propose  a  sequence  of  pulses  intended  to  preserve  the  state  of  a  qubit  in  the  presence  of  strong,  coherent 
coupling  to  another  quantum  system.  The  sequence  can  be  understood  as  a  generalized  SWAP  sequence  and  works 
in  formal  analogy  to  the  well-known  spin  echo.  Since  the  resulting  effective  decoherence  rate  of  the  qubit  state 
is  strongly  influenced  by  the  additional  system,  this  sequence  might  serve  to  protect  its  quantum  state  as  well  as 
negating  the  effects  of  the  coherent  coupling.  A  possible  area  of  application  is  in  large-scale  quantum  computing 
architectures,  where  spectral  crowding  of  the  resources  might  necessitate  a  method  to  mitigate  residual  couplings. 
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I.  INTRODUCTION 

Any  future  large-scale  quantum  computation  architecture 
will  involve  a  large  number  of  quantum-coherent  elements. 
Any  building  blocks  which  employ  similar  fundamental 
concepts  (i.e.,  superconducting  circuits,  quantum  dots,  ions, 
etc.)  will  necessarily  have  level  energies  similar  in  frequency 
to  many  others.  This  clustering  of  energies  may  lead  to 
unwanted,  spurious  couplings  between  different  elements, 
possibly  disturbing  the  operation  of  the  whole.  Proposed  ways 
around  this  problem  involve  tunable  coupling  elements  with 
large  on-off  ratios  [1,2]  as  well  as  hybrid  architectures  [3-5], 
where  elements  of  different  origins  and  thus  different  energy 
scales  are  used  for  different  tasks. 

In  this  article  we  propose  a  method  to  dynamically  negate 
the  effects  of  coherent  coupling  of  a  qubit  to  another  quantum 
system.  Our  proposal  consists  of  a  series  of  pulses  applied 
to  the  composite  system  and  is  inspired  by  the  well  known 
spin-echo  sequence  [6],  which  is  designed  to  mitigate  slow- 
noise  fluctuations  in  the  dynamics  of  a  qubit.  Since  in  the 
course  of  the  time  evolution  the  state  of  the  qubit  will  be  partly 
transferred  into  the  additional  quantum  system,  its  properties 
will  strongly  influence  the  resulting  state.  Specifically,  when 
considering  a  qubit  coupled  to  a  quantum  memory  element, 
we  find  an  improvement  in  the  relaxation  rate  of  the  qubit.  The 
sequence  might  therefore  also  serve  to  protect  the  state  of  a 
qubit  from  energy  relaxation. 

Schemes  working  towards  protecting  a  qubit  state  with  the 
use  of  quantum  memory  traditionally  focus  on  a  full  SWAP 
sequence  of  the  quantum  information  from  qubit  to  storage 
element.  The  information  then  resides  in  the  memory  for  the 
storage  time  T,  after  which  another  SWAP  sequence  is  used 
to  write  it  back  into  the  qubit;  see,  e.g..  Refs.  [3,7,8].  In 
contrast  to  this  traditional  approach,  the  sequence  we  propose 
only  partially  transfers  the  state  of  the  qubit  into  the  memory. 
The  possible  protection  is  therefore  weaker  than  in  the  usual 
approach,  but  on  the  other  hand  application  of  the  necessary 
pulses  can  be  significantly  faster  than  for  a  full  SWAP  sequence. 

Our  proposal  is  inspired  by  the  well-known  spin-  or  Hahn- 
echo  sequence  [6],  a  schematic  depiction  of  which  is  shown 
in  Fig.  1(a).  The  spin-echo  sequence  is  designed  to  protect 


a  two-level  system  (TLS)  in  a  Ramsey  experiment  from  the 
effect  of  slow,  environmentally  induced  fluctuations  of  its  level 
splitting.  It  finds  wide  application  in  NMR  and  more  recently 
the  field  of  quantum  computation  [9].  In  a  typical  Ramsey 
experiment,  a  slightly  detuned  7r/2  pulse  is  applied  to  a  two- 
level  system.  After  some  free  evolution  time  T  the  same  pulse 
is  applied  again  and  the  population  of  the  system  is  measured. 
The  signal  as  a  function  of  time  will  show  the  beating  of  the 
state  vector  against  the  rotating  frame  as  defined  by  the  pulse,  at 
a  frequency  given  by  the  pulse  detuning  56.  If  a  low-frequency 
noise  source  acts  on  the  value  of  the  level  splitting  6  of  the  two- 
level  system,  this  will  similarly  affect  the  detuning  56,  and  thus 
subsequent  runs  of  the  experiment  will  lead  to  different  results. 
This  phenomenon  is  known  as  pure  dephasing.  The  spin-echo 
sequence  now  consists  of  the  application  of  an  additional  tt 
pulse,  the  echo  pulse,  after  half  the  evolution  time  T /2.  Due 
to  the  effect  of  this  pulse,  the  state  vector  will,  during  the 
second  half  of  the  time  evolution,  effectively  retrace  the  path 
it  took  in  the  first  half.  The  system  will  at  time  T  thus  always 
reach  its  initial  state  again  and  the  effects  of  fluctuations  at 
frequencies  which  are  slower  than  the  inverse  evolution  time 
l/T  are  canceled  by  the  sequence. 

For  what  we  call  the  Ti-echo  sequence  we  want  to  apply  the 
same  principles  to  the  case  of  two  coupled  quantum  systems, 
namely,  a  qubit  and  a  quantum  memory.  In  the  case  of  strong 
coupling  between  the  two,  and  for  an  initial  state  where 
the  qubit  is  excited  while  the  memory  rests  in  its  ground 
state,  the  system  will  undergo  coherent  oscillations  where 
the  excitation  is  shifted  periodically  between  the  two  parts. 
The  frequency  of  these  oscillations  depends  on  the  strength 
of  the  coherent  coupling  as  well  as  on  the  detuning  between 
the  two  systems.  Assuming  both  these  quantities  to  be  constant 
in  time,  application  of  an  echo  pulse  after  half  the  evolution 
time  r/2  will  cancel  the  effect  of  the  coherent  oscillations. 
Since  we  have  a  coupled  system  to  which  we  apply  this  pulse 
sequence,  the  resulting  effective  decay  of  the  qubit  state  will 
be  strongly  influenced  by  the  decay  of  the  memory  [10].  In 
the  case  where  the  additional  quantum  system  shows  better 
coherence  properties  than  the  qubit,  this  will  allow  us  to 
find  an  effective  protection  of  the  qubit  state  from  energy 
relaxation. 
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FIG.  1.  Schematic  representation  of  the  pulse  sequence  for 
(a)  standard  spin-echo  sequence  and  (b)  the  proposed  Ti-echo 
sequence.  In  both  cases,  the  state  of  the  qubit  is  read  out  after  time  T. 
Also  in  both  cases,  the  application  of  echo  and  recovery  pulses  leads 
to  a  final  state  which  is  independent  of  T .  For  details,  see  the  text. 


Figure  1  (b)  gives  a  depiction  of  the  proposed  pulse  sequence 
and  compares  it  to  the  case  of  a  spin-echo  protocol.  In  the 
Ti-echo  sequence,  the  system  evolves  freely  for  a  time  r/2. 
After  this  period  of  free  evolution,  an  echo  pulse  applied, 
which  corresponds  to  a  tt  rotation  around  an  axis  perpendicular 
to  the  original  axis  of  rotation.  After  another  period  r/2  of 
free  evolution,  the  state  reached  by  the  system  will  no  longer 
depend  on  the  total  time  T  it  took  to  reach  it.  This  state  will  in 
general,  however,  not  be  the  starting  state  and  thus  has  to  be 
corrected  by  the  application  of  a  recovery  pulse  U r  .  The  exact 
nature  and  possible  realizations  of  the  pulses  and  Ur  will 
be  explained  in  detail  below. 


The  restriction  on  a  two-level  system  as  the  memory 
used  is  not  necessary  for  the  implementation  of  the  Ti-echo 
protocol;  we  use  Eq.  (1)  as  the  simplest  conceivable  model.  The 
original  motivation  for  these  work  came  from  research  with 
superconducting  qubits  coupled  to  intrinsic  two-level  defects 
[7,11,12].  This  system  naturally  realizes  the  Hamiltonian 
Eq.  (1)  and,  due  to  the  defects’  possible  long  coherence 
times  [11],  might  serve  as  a  testing  ground  for  our  proposal. 
Another  possibility  for  use  as  a  memory  would  be  a  harmonic 
oscillator  in  the  quantum  regime,  as  was  realized,  e.g.,  in 
Refs.  [12,13].  As  long  as  no  more  than  one  excitation  (i.e., 
either  one  photon  in  the  resonator  or  the  qubit  in  its  excited 
state)  in  total  is  involved  in  the  dynamics,  the  underlying 
Jaynes-Cummings  Hamiltonian  will  lead  to  the  same  results 
as  the  ones  presented  here.  Recent  experimental  efforts  have 
succeeded  in  coherently  transferring  quantum  information 
from  a  qubit  to  an  ensemble  of  spins  in  the  solid  state  [3,8,14], 
giving  yet  another  possible  experimental  realization  of  the 
model  Eq.  (1).  In  these  ensemble  systems,  the  weak  individual 
couplings  g  of  the  single  spins  act  collectively  to  realize  a 
strong  effective  coupling  proportional  to  \fN g  to  a  collective 
state  of  the  ensemble.  The  effective  Hamiltonian  is  again  of  the 
Jaynes-Cummings  type  and  thus  can  be  reduced  to  Eq.  (1)  for  a 
single  excitation.  We  want  to  note  that  there  are  other  transfer 
schemes  proposed  in  the  literature  which,  rather  than  using 
a  direct  existing  interaction  between  two  physical  systems, 
induce  such  an  interaction  by  use  of,  for  instance,  laser  or 
microwave  excitation;  see,  e.g..  Ref.  [15].  Eor  such  situations, 
our  scheme  seems  impractical  to  realize  since  the  interaction 
term  in  Eq.  (1)  will  be  necessary  at  all  times  for  our  proposal. 


II.  SYSTEM 

We  focus  our  theoretical  analysis  on  a  system  consisting 
of  two  coherently  coupled  two-level  systems;  one  is  the  qubit 
and  the  other  one  we  will  call  a  memory  element  from  now 
on.  The  Hamiltonian  reads 

Hd  =  -\iqOz  -  +  \v±ia+r_  +  ct_t+),  (1) 

where  a  and  r  are  the  Pauli  matrices  describing  qubit  and 
memory,  respectively.  The  qubit  has  a  level  splitting  of  size 
€q  while  the  memory  is  at  energy  6^,  and  they  are  coupled 
with  a  coupling  strength  i;^.  We  write  the  coupling  in  the 
rotating-wave  approximation,  assuming  6^,6^  ^  i;^.  We  also 
consider  the  regime  where  qubit  and  memory  are  close  to 
resonant,  €q  ~  6^ ,  meaning  that  the  coupling  between  the  two 
is  effective,  i.e.,  important  for  the  time  evolution. 

In  the  validity  range  of  the  rotating-wave  approximation,  the 
above  Hamiltonian  preserves  the  total  number  of  excitations. 
This  property,  when  decoherence  is  neglected,  enables  us  to 
effectively  decouple  the  dynamics  for  different  numbers  of 
excitation.  Eor  the  application  of  our  sequence,  we  assume  the 
memory  to  be  initially  in  its  ground  state  |0^ ) .  In  this  case,  and 
under  the  action  of  the  Hamiltonian  Eq.  (1),  in  order  to  show 
protection  of  the  qubit  state  it  is  sufficient  to  show  that  the  total 
ground  state  of  the  full  system  |0^,0^)  is  unaffected  and  that 
we  can  bring  the  state  back  to  itself.  We  will  not  be 

concerned  with  dynamic  or  geometric  phases  acquired  during 
time  evolution  since  they  can  be  negated  using  single-qubit 
gates. 


III.  EVOLUTION  ON  THE  EFFECTIVE  BLOCH  SPHERE 

To  illustrate  the  workings  of  the  Ti-echo  sequence,  we  will 
first  focus  on  the  case  of  purely  unitary  evolution,  i.e.,  zero 
decoherence.  In  this  case,  and  starting  with  the  qubit  in  its 
excited  state  and  the  memory  in  its  ground  state,  1 1^,0^),  the 
time  evolution  is  restricted  to  a  subspace  of  the  full  Hilbert 
space  containing  only  one  excitation.  We  can  describe  this 
subspace  by  an  effective  Bloch  sphere,  spanned  by  the  states 
containing  only  a  single  excitation  in  either  the  qubit  or  the 
memory,  i.e.,  span{|l^,0^),|0^,l^)}.  We  will  first  give  the 
expressions  for  U_i  and  U r  and  explain  their  effects  on  the  time 
evolution.  Later  in  this  section  we  will  also  show  how  to  realize 
these  pulses  without  the  need  for  high-frequency  controls, 
using  only  the  detuned  Hamiltonian  evolution  governed  by 
Eq.  (1). 

The  states  describing  the  effective  Bloch  sphere  are 


k+)  =  |i„o„)  ^  (0,0,1),  k-)  =  |o„u)  ^  (0,0,  -1), 


1 

7!^ 

1 

7!' 

1 

7!^ 

1 

7!' 


|x+)  =  ^(|1„0„>  +  |0„1„))  ^  (1,0,0), 


|x-)  =  — (|1^,0„)  -  |0,,1„))  ^  (-1,0,0),  (2) 


|y+>  =  —(|l„0„>+i|0„l„})^  (0,1,0), 


k-)  =  — (|l^,0„)-i|0^,l„))^(0,  -1,0), 
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where  the  indices  q  and  m  denote  states  of  the  qubit  and  the 
memory,  respectively,  and  the  vectors  (x,y,z)  describe  points 
on  the  pseudo-three-dimensional  Bloch  sphere. 

Restricting  the  Hilbert  space  to  the  one-excitation  subspace, 
the  time  evolution  operator  corresponding  to  Eq.  (1)  can  be 
written  as 


=  exp[-iHit} 


X 


/  -  cos  ^0 

V  sin^o 


(3) 


where  a>o  =^/vj_  +  Sco^,  tanfo  =  v±/Sa),  and  we  choose  for 
the  detuning  angle  ^  [0,^]-  Here,  Sco  =  €q  —  6^  is  the 
detuning  between  qubit  and  memory  and  we  write  Hi  for 
the  part  of  the  Hamiltonian  Hq  which  is  acting  on  the  one- 
excitation  subspace.  The  expression  Eq.  (3)  can  be  interpreted 
as  free  precession  on  the  Bloch  sphere  around  the  axis 
no  =  (sin §0,0,  cos  §o)  with  frequency  coq.  The  rotation  axis 
no  forms  an  angle  §o  with  the  z  axis.  Eor  zero  detuning, 
8(jo  =  0,  we  find  §o  =  7r/2  and  therefore  no  =  (1,0,0),  i.e., 
the  precession  takes  place  around  the  x  axis.  In  the  case 
of  strong  detuning,  8co  ^  v±,  we  have  §o  =  0  and  thus  the 
rotation  axis  is  no  =  (0,0,1).  If  the  initial  state  is  also  lying 
on  the  z  axis,  i.e.,  |1^,0^),  no  rotation  will  take  place.  This 
is  equivalent  to  saying  that  for  sufficiently  strong  detuning 
between  the  two  systems,  their  coupling  no  longer  affects  the 
time  evolution.  Einally,  for  nonzero  but  finite  detuning,  the 
resulting  rotational  axis  lies  somewhere  in  the  middle  between 
the  X  and  z  axes.  Eor  example,  for  8co  =  we  find  §o  =  7r/4 
and  no  =  ^(1,0,1).  Eor  our  choice  of  interaction  term  in  the 
Hamiltonian  Eq.  (1),  the  axis  around  which  rotations  take  place 
is  thus  fixed  on  the  xz  plane  of  the  effective  Bloch  sphere.  It  is 
important  to  note  that  this  is  no  special  property  of  our  choice 
of  interaction.  In  fact  any  physical  interaction  corresponding 
to  an  exchange  of  energy  between  the  qubit  and  the  memory 
will  give  a  similar  restriction  of  the  axis  of  rotation  to  a  plane 
in  the  effective-B loch- sphere  picture. 

In  keeping  with  the  ideas  of  spin  echo,  we  want  to  apply  a 
TT  pulse  around  an  axis  which  is  perpendicular  to  the  axis  no 
of  free  evolution,  and  we  want  to  apply  this  pulse  after  half  the 
total  evolution  time  T.  We  have  some  freedom  in  the  choice  of 
this  axis,  and  will  choose  it  lying  again  in  the  xz  plane  of  the 
one-excitation  Bloch  sphere.  This  selection  will  later  enable 
us  to  perform  this  rotation  using  only  detuned  Hamiltonian 
evolution,  without  the  need  for  high-frequency  controls.  The 
axis  n±  lying  on  the  xz  plane  and  perpendicular  to  no  is  simply 
nj_  =  (—cos §0,0,  sin§o).  An  operator  corresponding  to  a  tt 
pulse  around  n±  is  given  by 


UA^o)  = 


Vcc 


sin§o 

cos§o 


cos§o 

-sin§o 


(4) 


After  application  of  the  echo  pulse  U±,  another  period  of  free 
evolution  for  the  time  T jl  brings  the  system  into  the  state 
1 4^7),  which  is  independent  of  the  total  evolution  time  T .  The 
state  \^t)  before  the  application  of  the  recovery  pulse  Ur  is 


t  =  0  |10)  t  =  T/2  |10) 


EIG.  2.  (Color  online)  Step-by-step  illustration  of  the  Ti-echo 
sequence  on  the  effective  one-excitation  Bloch  sphere  for  the 
resonance  case  8a)  =  0,  with  a  starting  state  of  11^0^);  for  details 
see  the  text.  Red  arrows  show  the  path  of  the  state  vector  during  the 
two  time  spans  of  free  evolution;  the  blue  dashed  arrow  represents 
the  path  during  the  echo  pulse  U±.  The  recovery  pulse  is  not  shown, 
since  for  =  0  we  have  \^t)  =  —  i  I'ko)  -  The  illustration  depicts  a 
free-evolution  time  of  T  =  IttHScoq).  In  this  illustration  we  assume 
infinitely  sharp  pulses,  i.e.,  no  time  is  needed  for  the  application 
ofU^. 


given  by 

\^t)  =  Ui(T/2,^o)  Ui{T/2,^o)  |1?,0„) 

—  ^'(  sin§o|l^,0^)  “h  cos §0 10^ ,  )) 

=  (— sin2§o,0,  —  cos2§o),  (5) 

meaning  that  it  returns  again  to  the  xz  plane  of  the  effective 
Bloch  sphere.  It  is  important  to  note  again  that  the  state  1 )  is 
independent  of  the  evolution  time  T ,  i.e.,  the  echo  pulse  works 
as  expected. 

Eigure  2  shows  an  illustration  of  the  full  time  evolution 
without  application  of  the  final  recovery  pulse  U r  .  The  black 
arrow  in  the  pictures  shows  the  state  of  the  system  at  a  given 
time  t  when  starting  from  state  |1^,0^),  while  the  red  arrows 
depict  the  path  the  state  vector  took  on  the  Bloch  sphere  due  to 
free  Hamiltonian  evolution  in  the  initial  and  final  time  spans 
of  length  r/2.  The  blue  dashed  arrow  finally  depicts  the  path 
taken  during  application  of  the  echo  pulse  U±.  The  situation 
shown  in  Eig.  2  is  for  zero  detuning  between  qubit  and  memory, 
§0  =  7r/2,  where  \^t)  = 

To  recover  the  initial  state  \^o)  =  |l^,0m)  =  (0,0,1)  from 
\^t),  Eq.  (5),  there  are  again  several  possibilities.  Choosing 
again  a  rotation  axis  lying  on  the  xz  plane,  we  perform  another 
Tt  pulse  around  the  axis  n  r  lying  exactly  halfway  between  \  ^t) 
and  the  z  axis  on  the  xz  plane.  Eor  this  choice  we  find  riR  =  n± 
and  thus  Ur  =  U±. 

Eigure  3  gives  an  illustration  of  the  effect  of  the  full 
sequence  for  a  detuning  between  qubit  and  memory  of  = 
—  and  a  total  free-evolution  time  of  T  =  lit  jcoo.  The  time 
evolution  starts  and  ends  in  the  state  |1^,0^).  The  periods 
of  free  evolution  depicted  by  red  arrows  are  interspersed  by 
application  of  the  pulses  U±  and  Ur,  shown  in  blue. 
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EIG.  3.  (Color  online)  Illustration  of  a  full  Ti-echo  sequence  on 
the  effective  Bloch  sphere.  Red  arrows  show  the  path  of  the  state 
vector  during  the  free-evolution  periods,  while  the  dashed  blue  arrows 
show  the  effects  of  the  echo  pulse  and  the  recovery  pulse  Ur.  The 
black  vertical  arrow  depicts  the  state  of  the  system  at  the  start  and  end 
of  the  sequence.  The  depicted  situation  corresponds  io  8co  =  —Vj_, 
i.e.,  §0  =  371/4,  where  we  find  \^t)  =  k+)  and  a  total  evolution 
time  of  r  =  Itt/coq. 


The  pulses  U±  and  Ur,  given  by  Eq.  (4),  act  in  general 
on  both  qubit  and  memory  and  might  be  hard  to  realize 
experimentally.  However,  our  choice  of  rotation  axis  on  the 
xz  plane  for  both  pulses  enables  us  to  give  a  protocol  for 
the  Ti-echo  sequence  which  does  not  require  direct  unitary 
control  of  either  qubit  or  memory.  The  whole  sequence  can 
in  fact  be  realized  utilizing  only  Hamiltonian  evolution.  In 
this  scenario,  the  pulses  U±  and  Ur  correspond  to  periods 
during  which  the  qubit  is  tuned  to  a  different  working  point 
than  for  the  free  evolution  time  T.  The  working  point  for 
these  sequences  is  easy  to  calculate,  when  remembering  that 
the  pulse  we  want  to  effect  is  a  tt  pulse  around  an  axis 
which  is  perpendicular  to  the  original  free-evolution  axis  hq. 
The  position  of  no  is  completely  governed  by  the  detuning 
angle  =  Cretan  (u^ /5a;).  Therefore,  to  have  the  system 
evolve  around  an  axis  perpendicular  to  no,  we  simply  have 
to  change  the  detuning  8co  such  that  the  new  detuning  angle 


encloses  the  angle  it  j 2  with  §o-  We  find 

8(joi  =  —v\I8(d.  (6) 

In  order  to  generate  a  rotation  angle  of  it  we  have  to  leave 
the  qubit  at  this  detuning  for  a  time  U  =  it  jcoi,  where  coi  = 
is  the  oscillation  frequency  at  the  new  bias  point. 
For  these  parameters  of  detuning  and  evolution  time,  we  find 

ui(u,^i)  =  uA^oy 


Pu 


T  [units  of  1/71, g] 


P^ 


0  T/2  T  ^ 

EIG.  4.  (Color  online)  Illustration  of  the  Ti-echo  sequence  using 
only  detuned  Hamiltonian  evolution.  Both  pulses  and  Ur  can 
be  replaced  by  a  period  of  detuned  free  evolution  at  detuning  8co\  = 
—  v\/ 8 CO  for  the  time  U  •  The  red  line  depicts  the  qubit’s  level  splitting 
at  different  times.  The  resonance  condition  with  the  memory,  Cq  = 
is  indicated  by  the  dashed  line  in  the  middle.  The  situation  depicted 
corresponds  to  the  case  of  initially  8(0  =  where  8(Oi  =  —8(o. 


EIG.  5.  (Color  online)  Time  evolution  of  the  qubit’s  excited- state 
population  Pi^  as  a  function  of  free-evolution  time  T  with  (red  solid) 
and  without  (blue  dashed)  application  of  the  Ti-echo  sequence  and 
for  three  different  detunings  between  qubit  and  memory.  The  dotted 
green  line  shows  the  decay  of  the  qubit  itself  without  any  interactions 
with  other  quantum  systems.  The  pulses  are  assumed  to  be  infinitely 
sharp  in  time,  and  there  is  relaxation  acting  on  only  the  qubit,  not  the 
memory.  Parameters  are  (in  units  of  y\^q)  =  y^,^  =  0, 

=  5,  5a;  =  0  (top),  5  (middle),  15  (bottom). 
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Figure  4  shows  an  illustration  of  the  Ti  -echo  sequence  using 
only  Hamiltonian  evolution.  The  red  line  represents  the  bias 
point  of  the  qubit  as  a  function  of  time.  Resonance  with  the 
memory  is  indicated  by  the  dashed  black  line. 

We  want  to  stress  that  the  minimum  time  scale  for  the 
application  of  the  Ti-echo  sequence  is  given  by  2tjj;  =  In  I  coi. 
In  contrast,  the  minimum  time  scale  for  a  full  SWAP  operation 
(transferring  information  into  a  memory  and  back  again)  is 
24wap  =  27r/i;^,  which  is  always  longer  than  or  at  minimum 
the  same  length  as  2tj^ . 


IV.  EFFECTS  OF  DECOHERENCE 

Up  to  now  we  considered  the  ideal  case  of  purely  unitary 
evolution  governed  by  the  Hamiltonian  Eq.  (1).  In  a  real-world 
scenario,  however,  the  time  evolution  will  be  affected  not  only 
by  the  coupling  of  the  two  systems  to  each  other,  but  also  by 
their  inevitable  coupling  to  their  respective  environments.  To 
show  the  workings  of  the  Ti-echo  sequence,  we  therefore  have 
to  take  into  account  the  effects  of  decoherence  on  the  time 
evolution.  As  stated  initially,  the  sequence  has  the  properties 
that  it  might  serve  to  protect  the  qubit  state  against  relaxation 
in  the  specific  case  where  the  additional  coupled  system  is 
a  memory,  i.e.,  its  coherence  properties  are  better  than  those 
of  the  qubit.  We  will  therefore  focus  on  this  situation  in  all 
our  following  illustrations  although  our  theoretical  treatment 
is  valid  in  the  general  case. 

To  obtain  the  time  evolution  including  decoherence,  we 
solve  a  Bloch-Redfield  master  equation  [16,17]  in  secular 
approximation  [10]  for  the  density  matrix  p{t)  of  the  cou¬ 
pled  system  of  qubit  and  memory.  We  introduce  individual 
relaxation  rates  for  the  qubit  and  memory  as  Similarly 

we  write  the  pure  dephasing  rates  as  Yip^qjm  with  the  total 
dephasing  rate  given  by  72  =  Y^)• 

For  the  straightforward  case  where  the  pulses  JJ ^  and  JJr 
are  applied  as  simple  unitary  gates  with  no  time  cost,  we  are 
able  to  give  simple  and  intuitive  analytical  results  as  well  as  an 
analytical  solution  of  full  process  tomography,  demonstrating 
the  high  degree  of  protection  the  Ti-echo  sequence  offers 
against  energy  relaxation.  For  the  more  realistic  case  where  the 
pulses  are  effected  employing  purely  Hamiltonian  evolution, 
as  illustrated  in  Fig.  4,  we  will  present  numerical  results 
and  show  that  also  in  this  case  an  effective  protection  from 
decoherence  can  be  achieved. 
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For  clarity  we  will  first  focus  on  the  case  of  only  energy 
relaxation  acting  on  the  qubit  and  neglect  pure  dephasing. 
Then,  if  the  pulses  are  infinitely  strong,  i.e.,  it  does  not  take 
any  time  to  apply  them,  we  arrive  at  a  simple  solution  for  the 
density  matrix  p{t)  of  the  coupled  system.  For  a  starting  state  of 
1 ,0m)  we  calculate  the  population  of  the  qubit’s  excited  state 
=  ^[tr{p(r)(7^}  -h  1]  as  a  function  of  the  free-evolution 
time  T  and  get 

PuiT)  =  (7) 

where  y+  =  \{Yi,q  +  Ki,m)  is  the  mean  of  the  individual 
relaxation  rates  of  qubit  and  memory.  Notably,  this  result  is 
independent  of  the  detuning  Sco  between  qubit  and  memory. 

Figure  5  shows  the  time  evolution  of  the  qubit’s  excited- 
state  population  as  a  function  of  free-evolution  time  T 
from  analytical  calculations  when  considering  only  relaxation 
acting  on  both  subsystems.  The  blue  dashed  lines  show  free 
evolution  without  the  application  of  any  pulses.  The  red 
solid  lines  give  the  results  after  one  subjects  the  system  to 
the  Ti-echo  pulses.  With  the  sequence,  independent  of  5a;, 
a  simple  exponential  decay  is  seen,  following  Eq.  (7).  As 
a  comparison,  the  green  dotted  lines  shows  the  exponential 
decay  proportional  to  of  the  qubit  itself,  without  any 

coupling  to  another  quantum  system.  For  zero  detuning  (top 
panel),  the  free  evolution  without  any  pulses  shows  strong 
coherent  oscillations  of  the  excitation  between  qubit  and 
memory.  After  application  of  a  Ti-echo  pulse,  the  resulting 
decay  is  given  by  the  envelope  of  the  coherent  decay  curve 
(cf.  Ref.  [10]).  For  nonzero  detuning  (the  lower  two  panels), 
the  amplitude  of  the  oscillations  decreases  and  its  frequency 
(jOq  increases.  The  effective  decay  changes  from  decay  of  the 
mixed  systems  for  small  detuning  to  more  qubitlike  behavior 
for  stronger  detuning.  In  contrast,  the  decay  of  the  qubit  state 
after  application  of  the  Ti-echo  sequence  stays  the  same  and 
is  given  always  by  Eq.  (7). 

The  whole  sequence  including  decoherence  can  be  thought 
of  as  a  quantum  dynamical  map  s{pq)  acting  on  the  density 
matrix  pq  of  the  qubit.  We  can  write  this  map  in  terms  of  a 
basis  set  {£„}  of  all  operators  acting  on  the  Hilbert  space  of 
the  qubit  as 

~  ^  ^  Xmn  Em  Pq  (8) 

mn 

with  the  coefficient  matrix  x-  Choosing  as  basis  set  = 
{ 1 ,  (Tjc ,  }  we  get  for  the  matrix  x 


X 


1 

4 


/  1  -h  a  —  2a  d/2)K+7 

0 

0 

\l  sm€t 


0 

1  - 

/(I  -  a-^+0 
0 


0 

-/(I  -a-^+0 
1  - 

0 


(9) 


where  the  oscillation  frequency  is  6  ~  ~  6^  and  the  sole 

decay  rate  appearing  is  7+.  Again,  this  result  is  valid  for  only 
relaxation  acting  on  both  subsystems.  In  this  case  the  only 


relevant  decay  rate  in  the  time  evolution  is  the  mean  of  the 
individual  decay  rates  of  the  two  subsystems,  as  can  also  be 
seen  in  Fig.  5. 
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The  results  presented  up  to  now  are  for  the  ideal  case 
when  perfect  pulses  are  available  and  no  pure  dephasing  acts 
on  either  one  of  the  subsystems.  Above  we  have  given  the 
possibility  of  generating  the  Ti-echo  sequence  using  only 
detuned  free  evolution  without  the  need  for  high-frequency 
controls.  In  this  proposal,  however,  the  pulses  will  take  a 
finite  amount  of  time  during  which  decoherence  will  act 
on  both  systems  and  therefore  the  results  will  be  worse  than 
for  the  case  of  infinitely  short  pulses.  Additionally,  the  above 


pu 


t  [units  of  1/71, g] 


EIG.  6.  (Color  online)  Time  evolution  of  the  qubit’s  excited-state 
population  as  a  function  of  time  with  (solid  red)  and  without 
(dashed  blue)  application  of  the  Ti-echo  sequence  and  for  three 
different  detunings  between  qubit  and  memory.  The  pulses  are 
realized  in  this  case  using  a  detuned  Hamiltonian  evolution;  see  the 
text.  The  green  dotted  line  shows  again  simple  exponential  decay 
with  the  rate  Yi,q-  The  red  line  does  not  start  at  the  origin  due  to 
the  finite  pulse  time  2u  needed  for  application  of  the  sequence. 
Parameters  are  (in  units  of  y\^q)  y^p^q  =  0.5,  y\^rn  =  y(p,m  =  =  5, 

=  0  (top),  1  (middle),  2  (bottom). 


calculations  did  not  include  pure  dephasing  as  a  source  of 
decoherence. 

Figure  6  shows  the  time  evolution  of  the  qubit’s  excited- 
state  population  Pi^q  as  a  function  of  total  time  t  and  for 
different  detunings  when  the  sequence  pulses  are  realized 
by  a  detuned  Hamiltonian  evolution.  These  results  are  from 
numerical  calculations  solving  a  Lindblad  master  equation 
for  the  time  evolution  of  the  system’s  density  matrix.  Both 
relaxation  rates  yi  as  well  as  pure  dephasing  rates  have 
been  included  in  the  calculations.  The  red  solid  line  again 
shows  the  time  evolution  of  the  qubit’s  excited-state  population 
Pi^q  after  application  of  the  Ti-echo  sequence,  while  the 
dashed  blue  line  shows  the  free  evolution  without  any  pulses. 
The  green  dotted  line  is  again  for  comparison  with  simple 
qubit  decay  proportional  to  In  this  situation  the  qubit 

state  population  after  applying  the  pulses  still  shows  some 
residual  oscillations,  due  to  the  mixing  of  qubit  and  memory 
states  during  the  pulses.  The  overall  decay  is  somewhat  faster 
than  predicted  by  Eq.  (7),  mainly  due  to  the  additional  pure 
dephasing  acting  on  the  system.  With  the  exception  of  an  initial 
residual  oscillation,  which  is  suppressed  on  a  time  scale  given 
by  the  dephasing  rates  72,  the  probability  of  finding  the  qubit  in 
its  excited  state  Pi^q  is  a  smooth  function  of  time,  highlighting 
the  performance  of  the  Ti-echo  sequence.  Additionally,  in  the 
illustrated  case  of  interaction  with  a  memory  element,  the 
lifetime  of  the  qubit  state  is  enhanced  as  compared  to  the  free 
case. 

To  quantify  the  loss  of  qubit  coherence  during  the  appli¬ 
cation  of  the  echo  pulses.  Fig.  7  shows  a  plot  of  the  initial 
probability  of  finding  the  qubit  in  the  excited  state,  Pi,^(0), 


Fi,,(0) 


EIG.  7.  (Color  online)  Initial  decay  of  the  qubit  state  vector  due 
to  the  finite  pulse  length  obtained  from  numerical  calculations. 
The  red  solid  curve  shows  the  probability  of  finding  the  qubit  in  its 
excited  state  as  a  function  of  detuning  8(o,  after  only  the  pulses  have 
been  applied  to  the  system  via  the  detuned  Hamiltonian  evolution. 
The  total  time  that  has  passed  is  thus  2u,  while  no  free  evolution 
has  taken  place,  so  the  free-evolution  time  T  is  zero.  The  green 
dotted  curve  for  comparison  shows  the  value  of  ,  meaning  the 

value  of  qubit  excitation  probability  after  application  of  the  pulses 
without  any  influence  of  the  TLS.  The  blue  dashed  line  Anally  shows 
decay  with  the  composite  rate  y+.  One  can  see  that  for  small  initial 
detuning,  the  pulses  lead  to  decay  dominated  by  the  qubit  rate,  while 
for  stronger  initial  detuning  it  crosses  over  to  decay  of  the  coupled 
systems.  Parameter  values  are  (in  units  of  yi^q)  y^p^q  =0.5,  yi^m  = 
y<p,m  =  0,  =  5. 


032335-6 


Ti-ECHO  SEQUENCE:  PROTECTING  THE  STATE  OE  . . . 

when  no  free-evolution  time  T  has  passed,  i.e.,  only  the  pulses 
have  been  applied  via  detuned  evolution.  The  solid  red  line 
shows  Pi,^(0)  as  a  function  of  initial  detuning  Sco  including 
both  relaxation  and  dephasing  on  the  qubit.  For  comparison, 
the  dotted  green  line  gives  the  excitation  probability  of  the 
uncoupled  qubit  after  the  time  Itjj;  it  takes  to  perform  the 
gates.  Finally,  the  dashed  blue  line  gives  the  value  of  , 

i.e.,  the  effective  decay  of  the  resonantly  coupled  system.  As 
one  can  see,  close  to  resonance,  where  for  zero  time  T  the 
excitation  resides  mostly  in  the  qubit  during  the  pulses,  the 
decay  during  detuned  evolution  is  qubitlike,  i.e.,  proportional 
to  .  For  stronger  initial  detuning,  however,  the  detuning 

during  application  of  the  gates,  is  close  to  resonant 
and  therefore  the  effective  decay  becomes  similar  to  the 
decay  of  the  resonantly  coupled  systems  when  dephasing  is 
included  [10]. 

V.  DISCUSSION 

It  is  instructive  to  consider  the  limits  of  the  proposed  pulses. 
In  one  limit,  when  qubit  and  memory  are  initially  in  resonance, 
5a;  =  0,  the  pulses  correspond  to  a  very  strong  detuning  for  an 
equally  short  time.  This  effectively  leads  to  a  phase  gate  acting 
on  the  system,  which  in  the  effective-B loch- sphere  picture 
transfers  the  state  from  one  hemisphere  to  the  opposite  one  (cf. 
Fig.  2).  In  the  other  limit,  strong  initial  detuning,  5a;  ^  i;^, 
we  find  5a;i  =  0.  In  this  case  the  pulses  and  Ur  both 
correspond  to  a  full  transfer  of  quantum  information  between 
the  qubit  and  the  memory,  and  thus  to  a  full  SWAP  sequence 
between  qubit  and  memory.  This  situation  provides  an  intuitive 
explanation  for  the  resulting  effective  relaxation  rate  }/+  in 
Eq.  (7):  the  excitation  first  rests  in  the  qubit  for  a  time  T /2 
before  being  transferred  into  the  memory  and  resting  there 
for  another  time  span  T jl.  The  effective  decay  is  therefore 
governed  by  the  average  of  the  two  individual  rates.  In  the 
case  of  strong  detuning,  therefore,  use  of  a  traditional  SWAP 
sequence,  where  the  information  resides  for  the  full  time  T 
in  the  memory,  is  more  advantageous.  In  the  opposite  case, 
however,  the  Ti-echo  sequence  not  only  negates  the  effects 
of  coherent  coupling  but  also  provides  a  very  fast  way  of 
refocusing  the  qubit  state  and  offering  additional  protection 
from  relaxation. 

The  treatment  of  dephasing  above  was  done  by  introducing 
an  additional  pure  dephasing  rates  y^p  into  the  Redfield  tensor. 
A  more  realistic  treatment  has  to  take  into  account  the  origin 
of  pure  dephasing,  that  is,  slow  fluctuations  in  the  level  energy 
of  qubit  and  memory.  The  effect  of  these  fluctuations  on  the 
evolution  is  a  change  of  the  detuning  5a;  and  thus  of  the  angle 
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§0  as  well  as  the  oscillation  frequency  coq.  As  long  as  the 
evolution  time  T  is  shorter  than  the  inverse  frequency  of  the 
fluctuation,  the  Ti-echo  sequence  will  refocus  the  state  of 
the  qubit  also  in  this  case,  in  exact  analogy  to  the  Hahn 
echo.  For  the  special  case  of  zero  detuning  between  qubit 
and  memory,  5a;  =  0,  this  has  already  been  experimentally 
demonstrated  in  Ref.  [18].  In  that  work  it  was  shown  that 
additional  echo  pulses  at  equally  spaced  time  intervals 
can  provide  additional  protection  from  dephasing  [19].  This 
conclusion,  however,  is  true  only  for  the  slow-noise  spectrum 
leading  to  pure  dephasing.  Concatenation  of  pulses  will  then 
lead  to  a  further  decreased  sensitivity  of  the  system  to  the  low- 
frequency  parts  of  the  spectrum.  Considering  energy  relaxation 
into  a  Markovian  environment,  on  the  other  hand,  Eq.  (7) 
describes  the  maximal  achievable  protection  with  our  scheme. 

The  realization  of  the  Ti-echo  sequence  on  which  we 
have  focused  our  analysis  here,  using  detuned  Hamiltonian 
evolution,  is  only  one  of  many  possibilities.  In  concrete 
schemes  more  specific  sequences  might  be  more  practical. 
Depending  on  the  implementation  of  qubit  and  memory, 
two  qubit  gates  might  be  readily  available,  and  depending 
on  specific  gate  times  might  outperform  the  implementation 
without  high-frequency  controls  we  showed  here. 

VI.  CONCLUSIONS 

We  presented  the  Ti-echo  sequence  as  a  method  to 
dynamically  protect  the  state  of  a  qubit  from  the  effects  of 
interaction  with  another  coherent  quantum  system.  In  the  case 
where  the  additional  system  shows  better  coherence  properties 
than  the  qubit  itself,  it  additionally  serves  to  improve  the  qubit 
relaxation  time.  The  sequence  is  formally  similar  to  the  well- 
known  spin-echo  sequence  and  can  be  experimentally  realized 
without  the  need  for  high-frequency  controls.  Implementations 
in  different  systems  are  possible;  the  only  prerequisite  for  the 
proposed  implementation  is  the  tunability  of  the  qubit  tran¬ 
sition  frequency.  Possible  realizations  are  in  superconducting 
qubit  systems,  where  either  intrinsic  two-level  defects  [7]  or 
superconducting  resonators  [12]  would  play  the  role  of  the 
additional  elements,  and  implementation  is  straightforward. 

ACKNOWLEDGMENTS 

We  gratefully  acknowledge  discussion  with  D.  Pappas. 
M.  Sandberg,  and  M.  Vissers  and  want  to  thank  A.  Blais, 
G.  Grabovskij,  and  J.  Lisenfeld  for  their  comments  on  the 
manuscript.  We  acknowledge  financial  support  by  CFN  of 
DEG  as  well  as  NSERC. 


[1]  R.  A.  Pinto,  A.  N.  Korotkov,  M.  R.  Geller,  V  S.  Shumeiko,  and 
J.  M.  Martinis,  Phys.  Rev.  B  82,  104522  (2010). 

[2]  R.  C.  Bialczak,  M.  Ansmann,  M.  Hofheinz,  M.  Lenander, 
E.  Lucero,  M.  Neeley,  A.  D.  O’Connell,  D.  Sank,  H. 
Wang,  M.  P.  Weides  et  al,  Phys.  Rev.  Lett.  106,  060501 
(2011). 

[3]  Y.  Kubo,  C.  Grezes,  A.  Dewes,  T.  Umeda,  J.  Isoya,  H.  Sumiya, 
N.  Morishita,  H.  Abe,  S.  Onoda,  T.  Ohshima  et  at.,  Phys.  Rev. 
Lett.  107,  220501  (2011). 


[4]  Y  Kubo,  1.  Diniz,  A.  Dewes,  V  Jacques,  A.  Dreau,  J.  E.  Roch, 
A.  Auffeves,  D.  Vion,  D.  Esteve,  and  P.  Bertet,  Phys.  Rev.  A  85, 
012333  (2012). 

[5]  T.  Prey,  P.  J.  Leek,  M.  Beck,  A.  Blais,  T.  Ihn,  K.  Ensslin,  and 
A.  Wallraff,  Phys.  Rev.  Lett.  108,  046807  (2012). 

[6]  E.  Hahn,  Phys.  Rev.  80,  580  (1950). 

[7]  M.  Neeley,  M.  Ansmann,  R.  C.  Bialczak,  M.  Hofheinz,  N.  Katz, 
E.  Lucero,  A.  D.  O’Connell,  H.  Wang,  A.  N.  Cleland,  and  J.  M. 
Martinis,  Nat.  Phys.  4,  523  (2008). 


032335-7 


MULLER,  SHNIRMAN,  AND  WEIDES 


PHYSICAL  REVIEW  A  86,  032335  (2012) 


[8]  X.  Zhu,  S.  Saito,  A.  Kemp,  K.  Kakuyanagi,  S.-i.  Karimoto, 

H.  Nakano,  W.  J.  Munro,  Y.  Tokura,  M.  S.  Everitt, 
K.  Nemoto  et  al..  Nature  (London)  478,  221 

(2011). 

[9]  G.  Ithier,  E.  Collin,  R  Joyez,  R  J.  Meeson,  D.  Vion,  D.  Esteve, 
E.  Chiarello,  A.  Shnirman,  Y  Makhlin,  J.  Schriefl  et  al.,  Phys. 
Rev.  B  72,  134519  (2005). 

[10]  C.  Miiller,  A.  Shnirman,  and  Y.  Makhlin,  Phys.  Rev.  B  80, 
134517  (2009). 

[11]  J.  Lisenfeld,  C.  Muller,  J.  H.  Cole,  P.  Bushev,  A.  Lukashenko, 
A.  Shnirman,  and  A.  V.  Ustinov,  Phys.  Rev.  Lett.  105,  230504 
(2010). 

[12]  M.  Mariantoni,  H.  Wang,  T.  Yamamoto,  M.  Neeley, 

R.  C.  Bialczak,  Y.  Chen,  M.  Lenander,  E.  Lucero, 
A.  O’Connell,  D.  Sank  et  al..  Science  334,  61 

(2011). 


[13]  M.  Hofheinz,  E.  M.  Weig,  M.  Ansmann,  R.  C.  Bialczak,  E. 
Lucero,  M.  Neeley,  A.  D.  O’Connell,  H.  Wang,  J.  M.  Martinis, 
and  A.  N.  Cleland,  Nature  (London)  454,  310  (2008). 

[14]  J.  M.  Cai,  E.  Jelezko,  N.  Katz,  A.  Retzker,  and  M.  B.  Plenio, 
New  J.  Phys.  14,  093030  (2012). 

[15]  J.  J.  L.  Morton,  A.  M.  Tyryshkin,  R.  M.  Brown,  S.  Shankar, 
B.  W.  Lovett,  A.  Ardavan,  T.  Schenkel,  E.  E.  Haller,  J.  W.  Ager, 
and  S.  A.  Lyon,  Nature  (London)  455,  1085  (2008). 

[16]  E.  Bloch,  Phys.  Rev.  105,  1206  (1957). 

[17]  A.  G.  Redfield,  IBM  J.  Res.  Dev.  1,  19  (1957). 

[18]  S.  Gustavsson,  E.  Yan,  J.  Bylander,  F.  Yoshihara,  Y.  Nakamura, 
T.  P.  Orlando,  and  W.  D.  Oliver,  Phys.  Rev.  Lett.  109,  010502 
(2012). 

[19]  J.  Bylander,  S.  Gustavsson,  E.  Yan,  E.  Yoshihara,  K.  Harrabi, 
G.  Eitch,  D.  G.  Cory,  Y.  Nakamura,  J.-S.  Tsai,  and  W.  D.  Oliver, 
Nat.  Phys.  7,  565  (2011). 


032335-8 


